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ABSTRACT 

Aims. A high number of embedded clusters is found in the Galaxy. Depending on the formation scenario, most of them 
can evolve to unbounded groups that are dissolved within a few tens of Myr. A systematic study of young stellar 
clusters showing distinct characteristics provide interesting information on the evolutionary phases during the pre-main 
sequence. In order to identify and to understand these phases we performed a comparative study of 21 young stellar 
clusters. 

Methods. Near-infrared data from 2MASS were used to determine the structural and fundamental parameters based on 
surface stellar density maps, radial density profile, and colour-magnitude diagrams. The cluster members were selected 
according to the membership probability, which is based on the statistical comparison with the cluster proper motion. 
Additional members were selected on basis of a decontamination procedure that was adopted to distinguish field-stars 
found in the direction of the cluster area. 

Results. We obtained age and mass distributions by comparing pre-main sequence models with the position of cluster 
members in the colour-magnitude diagram. The mean age of our sample is ^ 5 Myr, where 57% of the objects is 
found in the 4-10 Myr range of age, while 43% is < 4 Myr old. Their low E(B-V) indicate that the members are not 
suffering high extinction (Ay < 1 niag), which means they are more likely young stellar groups than embedded clusters. 
Relations between structural and fundamental parameters were used to verify differences and similarities that could be 
found among the clusters. The parameters of most of the objects have the same trends or correlations. Comparisons 
with other young clusters show similar relations among mass, radius and density. Our sample tends to have larger 
radius and lower volumetric density, when compared to embedded clusters. These differences are compatible with the 
mean age of our sample, which we consider intermediary between the embedded and the exposed phases of the stellar 
clusters evolution. 
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1. Introduction 

It is generally known that most stars are formed in groups 
or clusters. However, detailed studies on the initial pro- 
cesses of s tar formation are restricted to isolated dense cores 
of clouds (|Shu et al. 19871 l200lj) . 

The first stages of multiple star formation is usually 
evaluated through millimetric surveys, which are successful 
probing the scenario of stellar clusters formation. Based 
on 1.2mm data of the Mon OBI region, for instance, 
|Peretto et al. (200*5) discovered 27 proto-stellar cores with 
diameters of about 0.04 pc and masses ranging from 20 to 
40 M associated to the young stellar cluster NGC 2264. 
These results reveal the physical conditions for multiple 
massive star formation in a region that shows a wide range 
of masses (Dahm 2008|) and ages (|Flaccomi o et al. 1 997, 
Reb ull et al. 2002)) indicating the occurrence of a large va- 
riety of young stellar groups. 

In addition to millimetric studies, the evolution dur- 
ing the pre-main sequence (pre-MS) is more closely sur- 
veyed by using infrared data. Particularly, near-infrared 
(NIR) provides information about the circumstellar struc- 
ture of the cluster members, whose physical conditions are 
direc tly related to the pre-MS ev olutionary stage of the 
star. i Gregorio-Hetem et al. (2009)| for instance, have used 



X-ray results of sources detected in CMa Rl, combined 
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with NIR and optical data, to identify the young popu- 
lation associated to this star-forming region. They studied 
two young clusters having similar mass function, but one 
of them is older than the other by at least a few Myr. A 
mixing of populations seems to occur in the inter-clusters 
region, possibly related to sequential star formation. 

Systematic studies of young stellar clusters can directly 
probe several fundamental astrophysical problems, such as 
formation and evolution of open clusters, and more general 
problems, like the origin and ear ly ev olution of sta rs and 
planetary systems ([Adams et al.l 2006, Adams 2010). 

Kharche nko et all (2005) presented angular sizes of clus- 
ter core and "coronae" for a sample of 520 Galactic open 
clusters, and Piskuno v et al.l (2007) revised these data to 
obtain tidal radius and core size estimated by fitting King 
profile to the observed density distribution. 

Carpenter (2000) used 2MASS data to derive surface 
density maps of stellar clusters associated to molecular 
clouds, determining cluster radius and number of members 
based on the distributed population. 

I Adams et al.l (2006) explored the relations between clus- 
ter radius and stellar distribution in order to study the 
dynamical evolution of young clusters. Numerical simu- 
lations were used to reproduce different initial conditions 
of cluster- format ion. The statistical calculations developed 
by I Adams et al.l (2006) were based on cluster membership 
and cluster radius correlations observed in the samples pre- 
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sented by Car penter] (2000) and Lada & Lada (2003), who 
compiled the first extensive catalogue of embedded clusters 
providing parameters such as mass, radius and number of 
members. 

Lada & Lada (2003), hereafter ILL03[ verified the occur- 
rence of a high number of embedded clusters within molec- 
ular clouds in the Galaxy. However, most of these objects 
probably lose their dynamical equilibrium that would dis- 
solve the group turning it into field- stars. 

120091 studied the evolution of a sample of 23 massive 
clusters younger than 20 Myr, suggesting two distinct se- 
quences. Depending on size and density, a bimodal distribu- 
tion is found for the cluster size as a function of age. Based 
on the mass-radius and density-radius dependence observed 
in the sample from ILL031 Pfalzner (2011) proposed a sce- 
nario of sequential formation of "leaky" (exposed) massive 
clusters. 

On the other hand, differently of the time sequence pro- 
posed by Pfalzner (|201ip , the relations of clust er properties 
can be interpreted as formation condition ( Ada ms et al.l 
2006. lAdams 2010p . 

The structure of embedded clusters traces the physical 
conditions of their star-forming processes, since the origin 
and the evolution of stellar clusters are related to the dis- 
tribution of dense molecular gas in their parental cloud. 
ILL03I proposed two basic structural types, according to the 
cluster surface density. Classical open clusters are charac- 
terised by a high concentration in their surface distribution, 
whose radial profile is smooth and can be reproduced by 
simple power law functions or King-like (isothermal) po- 
tential. This type of embedded cluster is considered cen- 
trally concentrated. On the other hand, hierarchical type 
clusters exhibit surface density distribution with multiple 
peaks and show significant structure over a large range of 
spatial scale. Although there are clear examples of both 
types of structures, their relative frequency is unknown. 

The main goal of the present work is to characterise 
a large sample of young stellar clusters and to perform a 
comparative study aiming to verify their similarities and 
differences, which are related to their evolutionary stages 
in the pre-MS. 

The methods adopted by us have been commonly used 
in t he characterizati on of open clusters like those developed 
by IBonatto &; Bical (2009a), for instance. Stellar density 
maps are built from NIR data in order to derive parame- 
ters on the basis of radial density profile. These parameters 
define the structure of the cluster, which is related to its 
origin and evolution. 

The analysis based on radial density profile is unprece- 
dented to 86% of our sample. Considering the lack of 
systematic studies comparing pre-MS stellar clusters, dis- 
tributed in different galactic regions, the present work aims 
to provide sets of structural and fundamental parameters, 
determined in a uniform data analysis that may contribute 
to the discussion on the origin of stellar groups. 

In Sect. 2 we describe the sample by presenting the se- 
lection criteria and the decontamination method for dis- 
tinguishing cluster members from field-stars. Section 3 is 
dedicated to determine the structural parameters, which 
are used to accurately determine distance, age and mass 
that are presented in Sect. 4. A comparative analysis is 
performed in Sect. 5, while Sect. 6 summarises the results 
and the main conclusions. Finally, Appendix A displays all 
the plots of the entire sample. 



Loktin (19 94) 
2005), IPiskunov et all 



2. Description of the sample 

Several compilations of open stellar clus ters are avail 
able in the literature, for examp le |Lynga ( 1987) 
Mermilliod (1995)|ILL0^liaiarchenko et at 
WEBDAQ and DAMlE 
Dias et al. 20021 |2006|) . among others. In Sect. 2.1 
we present the criteria used to select the clusters, while 
Sect. 2.2 describes the method used to exclude field- stars 
that were found in the same direction of the cluster. 



2.1. Clusters selection and observational data 

Aiming to focus our study on pre-MS objects, the DAML 
catalogue was used to select young stellar clusters with ages 
in the 1-20 Myr range. Distances smaller than 2 Kpc were 
used as selection criterion in order to ensure good quality of 
photometric data. Despite our selection is limited to south- 
ern objects, they are distributed in different star- forming 
regions enabling us to compare diverse environments. Table 
1 gives the list of 21 selected clusters. 

A membership probability (P) estimated on basis of 
proper motion is given by DAML. For the selection of stars 
belonging to the clusters, the values P> 50% were adopted 
to indicate the possible members, hereafter denoted by P50. 

In order to confirm the cluster centre coordinates avail- 
able in the literature, we evaluated the distribution of num- 
ber of stars as a function of right ascension (a) and declina- 
tion (5). Gaussian profiles were fitted to these distributions, 
in order to estimate their centroid. A good agreement with 
the literature was found, within the errors estimated by 
the fitting that are AS ~ 1 arcmin and 0.75'< Aa < 2.25'. 
Table 1 gives the error on a indicated in between parenthe- 
sis. 

The NIR data, which provide maximum variation 
among colour-m agnitude diagrams (CMDs) of c lusters with 
different ages (Bonatto, Bica & Girardi 2004), were ex- 
t racted from 2MAS S (All Sky Catalogue of Point Sources 
- ICutri et al. 20 03). The JHK S magnitudes were searched 
for all stars located within a radius of 30 arcmin from the 
cluster centre. Only good accuracy photometric data were 
used in our analysis, by selecting objects with AAA qual- 
ity flags that ensure the best photometric and astrometric 
qualities ([Lee et al. 20 05). 



2.2. Field-stars decontamination 

Most of our clusters are projected against dense fields of 
stars in the Galactic disk, making difficult to distinguish the 
cluster members. The first step to determine the structural 
and fundamental parameters is to proceed the field-stars 
decontamination. 

This process is based on statistical analysis by compar- 
ing the stellar density of the cluster area with a reference 
region, near to the cluster, according to their characteristics 
in the CMD. The decontamination algorithm was employed 
as follows: 

— the whole range of magnitudes and colours is divided 
into three dimensional cells (J, J-H and J-Ks); 



1 http://www.univie. ac. at/webda/ 

2 http://www. astro. iag. usp. br/^wilton/ 
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Table 1. List of clusters and their structural parameters. 
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(h m s) 
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(pc) 


(PC ) 
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(PC ) 
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Collinder 205 


09 00 31(5) 


-48 59 


3.4±1.2 


9±3 


51±16 


4.8±1.9 


0.30±0.03 


6.6±2.6 


0.09±0.03 


?? 


Hogg 10 


11 10 40(9) 


-60 14 


1.9±0.4 


16±2 


9±4 


7.6±1.7 


0.46±0.28 


1.6±0.3 


0.24±0.15 


H 


Hogg 22 


16 46 35(6) 


-47 05 


1.9±0.5 


20±4 


7±3 


9.0±3.8 


1.91±1.10 


1.3±0.2 


1.03±0.65 


H 


Lynga 14 


16 55 03(5) 


-45 14 


1.1±0.4 


55±18 


153±52 


16.7±7.0 


0.22±0.04 


3.8±1.3 


0.20±0.08 


CC 


Markarian 38 


18 15 15(5) 


-19 01 


1.6±0.5 


18±6 


44±15 


6.8±2.3 


0.13±0.03 


3.4±1.2 


0.08±0.03 


CC 


NGC 2302 


06 51 55(4) 


-07 05 


2.3±0.5 


8±2 


16±4 


4. Oil. 4 


0.54±0.10 


3.0±0.6 


0.23±0.06 


CC 


NGC 2362 a 


07 18 42(4) 


-24 28 


1.9±0.4 


13±3 


62±14 


3.3±1.0 


0.26±0.05 


5.8±1.7 


0.13±0.04 


CC 


NGC 2367 


07 20 06(4) 


-21 53 


2.9±0.7 


5±1 


13±3 


2.3±0.6 


0.34±0.11 


3.6±0.7 


0.12±0.05 


?? 


NGC 2645 


08 39 05(6) 


-46 14 


2.1±0.5 


13±3 


48±12 


7.5±2.1 


0.22±0.04 


4.8±1.2 


0.10±0.03 


?? 


NGC 2659 


08 42 37(6) 


-44 59 


3.6±0.5 


9±1 


11±2 


5.2±0.9 


1.93±0.37 


2.2±0.3 


0.53±0.12 


H 


NGC 3572 


11 10 26(8) 


-60 15 


2.5±0.5 


19±3 


76±15 


7.7±1.9 


0.16±0.03 


5.0±1.1 


0.06±0.02 


H 


NGC 3590 


11 12 59(7) 


-60 47 


1.5±0.4 


30±7 


96±27 


11.7±4.5 


0.29±0.06 


4.2±1.2 


0.20±0.07 


CC 


NGC 5606 


14 27 48(8) 


-59 38 


2.6±0.4 


11±2 


25±5 


4.8±0.9 


0.51±0.13 


3.2±0.6 


0.20±0.06 


CC 


NGC 6178 


16 35 47(6) 


-45 39 


1.7±0.4 


24±4 


32±11 


12.2±2.7 


0.30±0.12 


2.3±0.5 


0.18±0.08 


H 


NGC 6604 b 


18 18 04(3) 


-12 15 


2.6±0.7 


11±3 


17±5 


4.4±1.3 


0.50±0.11 


2.6±0.6 


0.19d=0.07 


H 


NGC 6613 


18 19 58(4) 


-17 06 


2.9±0.7 


9±2 


82±19 


4.9±1.2 


0.12±0.01 


10.6±3.5 


0.04±0.01 


CC 


Ruprecht 79 


09 40 59(7) 


-53 51 


3.9±0.5 


7±1 


4±1 


3.6±0.5 


2.2±0.8 


1.5±0.2 


0.57±0.22 


CC 


Stock 13 c 


11 13 04(8) 


-58 53 


2.3±0.3 


15±2 


42±14 


3.8±0.9 


0.11±0.05 


3.8±1.0 


0.05±0.02 


H 


Stock 16 


13 19 30(8) 


-62 38 


2.8±0.7 


15±3 


8±4 


5.8±1.6 


1.2±0.6 


1.5±0.3 


0.42±0.24 


H 


Trumpler 18 


11 11 27(9) 


-60 40 


4.6±0.8 


9±1 


15±4 


2.5±0.4 


0.38±0.10 


2.8±0.5 


0.08±0.03 


H 


Trumpler 28 


17 37 02(5) 


-32 38 


1.4±0.5 


9±4 


25±9 


12.3±4.5 


0.84±0.16 


3.7±1.4 


0.61±0.24 


H 



Columns description: (1) Identification; (2, 3) J2000 coordinates; (4) cluster radius; (5) background density; (6) core density; (7) observed average density; 
(8) core radius; (9) density-contrast parameter; (10) ratio of core-size to cluster-radius; (11) cluster type: hierarchical (H), centrally concentrated (CC), 
undefined (??). Notes: parameters from lPiskunov et al.l (a) i? c =0.6±0.4 pc; (b) R c = 2.4±1.4pc; (c) J R C = 1.0±0.8 pc. 



for each cell (i), the stellar density of field-stars (erf = 
iif I a/) is obtained by counting the number of stars 
(iif) appearing in the reference region (a/); 
similarly, the total stellar density (cr t = n^ / a c ) is ob- 
tained by counting all stars (n t ) in the cluster area (a c ); 
the number of field-stars (rif), appearing in the cluster 
area, that have colours and magnitude similar to those 
estimated in the reference region is calculate by: 



rif = ^7 xn t 
at 



(i) 



Finally, the number of field-stars in the cluster area 
is randomly subtracted in each cell, remaining TV = 
Si (n ti — rif.) members of the cluster. 



In order to minimize errors introduced by the param- 
eters choice and uncertainties on 2MASS data, the decon- 
tamination algorithm is applied to several grids by adopting 
three different positions in the CMD. By defining a cell size 
AJ = 0.5 mag, the algorithm uses a grid starting at Ji 
and two other grids starting at Ji ± |AJ. This dithering 
is also made for both colours, A(J-H) = A(J- K5) = 0.1 
mag, providing 27 different results for each star. By this 
way, we obtain the probability that a given object should 
be considered a field-star. All possible members (P50) have 
0% of probability of being a field-star. 

We considered as candidate members (denoted by P?) 
the stars that were not removed by the field-star decontam- 
ination. Both, P50 and P? members were used, but having 
different weight, in the estimation of distance and age, pre- 
sented in Sect. 4.2. 

The adopted decontamination procedure has been used 
in studies of several types of cluster s: objects showing low 
contrast relative to the field-stars (IBica &; Bonatto 20051 
embedded 



clusters (Bonatt o, Santos Jr. fc Bica 2006), 
young clusters ([Bonatto et al. 2006|h among others. 



3. Structural parameters 

The structural parameters were determined on basis of 
stellar surface density, derived from stellar surface-density 
maps and radial density profile that are detailed in the 
Sects. 3.1 and 3.2. 

The first step, in the calculation of the stellar surface 
density, is to enhance the contrast between surface distri- 
bution of cluster members and field-stars by using a colour- 
magnitude filter. 

Pre-MS isochrones were adopted to establish the colour- 
magnitude filter limiting a CMD region that should con- 
tain only cluster members. Once the observed magni- 
tudes were unreddened, using the visual extinction given 
in the literature, we disregarded all the objects lying out 
of the range denned by the colour-magnitude filter. This 
filter is successful to accentuate the structures and to re- 
duce the fluctuations caused by the presence of field-stars 
Pica & Bonatto 2005)) . 

3.1. Stellar surface-density maps 

We obtained spatial distribution maps of the stellar surface 
density (a) given by the number of stars per arcmin 2 for 
the clusters and their surroundings. Appendix A presents 
the results for the whole sample, while Fig. 1 shows the 
map derived for the cluster Lynga 14, as illustration. The 
left panel shows the entire studied area, whose surface den- 
sity was calculated in cells of \A(acos(S cente r))\ = |A5| = 
2.5 arcmin 2 , where Aa and AS are the steps on right as- 
cension and declination, respectively. The dashed annulus 
corresponds to the reference field-stars region and the cen- 
tral circle indicates the cluster area. 

A zoom view (Fig. 1 central panel) displays a more de- 
tailed structure that was obtained by using smaller cells 
(1.0 arcmin 2 ) for the star counting process. Following the 
classification proposed by ILL03| the surface density maps 
were visually inspected to characterise the clusters. Last 
column of Table 1 indicates if the type is Hierarchical (H) 
or Centrally Concentrated (CC), according to the definition 
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A[occos(5)] A[acos(5)] 1 R (arcn 



Fig. 1. Left: Stellar surface-density map (a (stars/arcmin 2 )) obtained for the region of 30 arcmin around Lynga 14. The 
comparison field-stars area is indicated by dashed lines, while the full line indicates the cluster area. Center: a zoom of 
the <j map indicating by crosses the position of objects with membership probability P>70%. Right: The distribution of 
the stellar density as a function of rad ius. The best fitting of observed radial density profile, indicated by the full line, 
was obtained by using the model from |King (1962)[ A red line indicates the background density (cr^). 



discussed in Sect. 1. About 38% of the sample (8/21) shows 
a single major peak of density being classified by CC. Half 
of the sample (10/21) has multiple peaks and were consid- 
ered of H type, while three objects presenting filamentary 
density distributions, have undefined type (marked as "??" 
in Table 1). 



3.2. Radial density profile 

Aiming to quantify the stellar distribution of the clusters, 
we evaluated their radial density profile (RDP) by using 
concentric rings to calculate the surface density. 

Some of the structural parameters were obtained by fit- 
ting a theoretical RDP to the observed data. The adopted 
function is similar to the empirical model from King (1962)| 
given by: 



o-(r) = cr bg + 



1 + (r/r c y 



(2) 



where <r(r) is the stellar surface density (stars/arcmin 2 ), r 
is the radius (arcmin) of each concentric annulus used in the 
star counting, ctq and r c are respectively the density and the 
radius of the cluster core. The average density measured 
in the reference region (dbg) was calculate separately, in 
order to diminish the number of free parameters in the RDP 
fitting. 

Figure 1 (right panel) shows an example of observed 
radial distribution of surface density and the best fitting of 
the King's profile, which was obtained by adopting the x 2 
method based on the parameters cr and r c . 

Aiming to determine the cluster radius (R) more accu- 
rately, we verified the point where the cluster stellar density 
reaches the background density. Figure 1 displays a red line 
corresponding to a = a^g that was used to find R. 



A quantitative estimation on how compact is the clus- 
ter can be obtained from the density-contrast parameter 
proposed bv lBonatto fc Bical (2009a): 



5 C = 1 + ^- 



(3) 



where compact clusters have 7 < S c < 23. According this 
this criterion, only NGC 6613 can be considered compact 
(S c ~ 10). 

Two other parameters are usual in comparative analy- 
sis: the average density (< n* >), calculated by dividing 
the total number of observed members by the cluster area, 
and the ratio of core-size to cluster-radius (r c /R). We could 
expect small values of r c /R for the younger clusters, since 
their members would not have had enough time to disperse 
away from the centre. 

Table 1 gives the structural parameters and the uncer- 
tainties derived from the RDP fitting. 

Three of our objects are foun d in the cat alogue of 
Galactic open clusters presented by Pis kunov et al.l Based 
on proper movement, they identified bright stars ( V < 
14 mag) and used the radial profile of the region con- 
taining possible members (P=14-61%) to define the "coro- 
nae" radius of the cluster, while the concentration of 
the probable mem bers (P>61%) defines the core radius 
(Kharche nko et all 2005). Therefore, their sample is 5 to 
10 times less numerous than our sample that includes low- 
mass stars detected by 2MASS. Besides, we focused on the 
main stellar distribution of the cluster (radius < 6'), while 
they studied larger areas (radius ~ 15'), seeking for the 
tidal radius (Rt)- 

These different definitions of cluster radii, imply in sys- 
tematically smaller values when comparing our results with 
the structural parameters listed by them. However, these 
results cannot be directly compared because we are not 
dealing with the same kind of stellar groups. 
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In fact, the tidal radius obtained by Pisk unov et al.l 
(2007) for NGC 6604 (i? t =8.8pc); NGC 2362 (6.2pc); and 
Stock 13 (7pc) is more compatible with the tidal radius 
that we estimated for these clusters: r t = 7.7±2.3; 8.1±2.4; 
and 7.0±2.1, respectively, by adopting the expression used 
by Saurin, Bica & Bonatto (|20T2j) : 



M 

Mgai 



L / 3 d, 



GC 



(4) 



where M is the cluster mass, M ga i is the Galactic mass and 
doc is the Galactocentric distance, given by: 



gal 



Vgc d GC 
G 



(5) 



where Ygc — 254±16 km/s is the circular rotation velocity 
of the Sun at R GC = 8.4±0.6 kpc (jReid et al. 2009j) . 

On the other hand, the R ex N 0,5 dependence dis- 
cussed in Sect. 5.1 distinguishes the sample studied by 
Pis kunov et al.l (2007) from other young clusters present- 
ing size and membership relations similar to our sample. 
It seems to be more consistent comparing the cluster ra- 
dius that we o btained with the core radius obtained by 
IPiskunov et al.| given as notes in Table 1. 
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Fig. 2. Left: Colour-colour diagram for Lynga 14. The MS 
and the ZAMS are indicated by full lines, while the locus 
of giant star s is represented by a dotte d line. Reddening 
vectors from Rieke and Lebofsky (1985)] are shown by dot- 
dashed lines. Right: Colour-magnitude diagram showing the 
isochrones and evolutionary pre-MS tracks from lSiess et al.l 
(2000). Cluster members are indicated by open diamonds 
(P50) and dots (P?). 



4. Revisiting the fundamental parameters 

Adopting the same procedure described in Sect. 2.2, the 
field-stars decontamination was refined by using the accu- 
rate estimative of cluster radius. The size of the reference 
region was also re-evaluated according to the cluster area, 
allowing us to more closely define the sample of candidate 
members. 

In Sect. 4.1 we adopt the extinction available in the 
literature to correct the observed colours, which were fitted 
to the MS intrinsic colours. An iterative fitting process was 
adopted to accurately determine E(B-V). 

Despite the distance and age of our clusters are available 
in the literature, these parameters were checked by us on 
the light of the well determined structural parameters, as 
described in Sect. 4.2. 



4.1. Colour excess 

The (J-H) and (H-K s ) colours were used to estimate the 
extinction, by evaluating the position of the cluster mem- 
bers compared to the MS. In Fig. 2 (left panel), we 
plot the intrinsic colour s of MS and giant stars given by 



Bessel and Brett (1988)| as well as the corresponding red- 
dening vectors given by |Rieke and Lebofsky (1985) For 
comparison , the Zero Age Main- Sequence (ZAMS) from 
Sies s et aLl (2000) is also plotted, which has quite the same 
distribution of the MS, mainly for massive stars. 

Adopting the normal extinction law Ay = 3.09E(B — V) 
from Savage & Mathis (1979) and the relation 4^- from 



Cardelli et al. (1989) , the observed colours of the cluster 
members were unreddened and fitted to the MS intrinsic 
colours. This fitting is based on massive stars, mainly P50 
members. Table 2 gives the E(B-V) that provides the best 
fitting, which is in good agreement with those available in 
the literature, within the estimated errors. Only a few ob- 
jects had E(B-V) incompatible with the literature. In these 



cases, the procedure described in Sect. 3 for field-stars de- 
contamination by using colour-magnitude filter was reap- 
plied with the extinction derived by us, and the structural 
parameters were determined for this refined sample of clus- 
ter members. 

The (J-H) and (H-K S ) G colours were also used to reveal 
the stars with K-band excess. Young clusters are expected 
to have a large number of members showing high E(H-K S ) 
([Lada et al. 1996[) . These stars appear in the right side of 
the MS reddening vector in the colour-colour diagram (Fig. 
2). The fraction fx was calculated by dividing the number 
of stars having large (H-K s ) by the total number of cluster 
members. 

4.2. Evaluation of distance and age 

The unredenned magnitude and colour J x (J-H) G of 
the cluster members were compared to pre-MS mod- 
els ([Siess et al.l 2 000) and also MS Padova models 
([Girard et al. 2 002) that were required to fit the colours 
of massive stars. 

Figure 2 (right panel) shows the CMD obtained for 
Lynga 14, for illustration. The distances were confirmed 
by searching for (m — M)j Q that best fits the position of 
massive stars. The error bars were estimated from the mini- 
mum and the maximum distance module that provide good 
MS fittings. Within the uncertainties, the distances esti- 
mated from the MS fitting are in good agreement with the 
literature (differences are lower than 30%), excepting for 
Trumpler 18. Our results are presented in Table 2 and were 
used to convert into parsec the angular measurements of 
the structural parameters. 

The number of stars as a function of age was obtained by 
counting the objects in between different pairs of isochrones 
in the CMD. We use bins of 5Myr, excepting the two first 
that correspond to 0.2 - 1 Myr and 1-5 Myr ranges. For 
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Fig. 3. Left: Age distribution of Lynga 14 members (thick 
line) showing the contribution of P50 (dotted line) and P? 
(dashed line) objects. Right: Observed mass distribution in- 
dicated by crosses with error bars. The thick line represents 
the mass function 0(m) fitting. 



MS objects (above the 7 M track) we adopted the age 
estimated by fitting the data with Padova model. Figure 
3 (left panel) shows the histogram of age distribution as a 
function of fractional number of members, where possible 
members (P50) and candidate members (P?) are displayed 
separately. Figures A3 to A7 show the plots used to evaluate 
age and mass for the whole sample. 

The determination of age on basis of these histograms 
may not be obvious. Some of the clusters have members 
separated in two different ranges, making difficult the choice 
between the two options. The first choice of age is based on 
the most prominent peak formed by the P50 members. If 
the second peak has a similar number of members, we adopt 
the smaller value of age that is listed in Table 2, while the 
second option is given as notes in the same table. 

Due to the large uncertainties on the age estimation, 
our conclusions are based on mean values. Considering that 
12/21 objects of the sample (57%) are in the range of 4 - 10 
Myr, and 9/21 (43%) are younger than 4 Myr, we suggest 
for our clusters a mean age of ~ 5 Myr. 

4.3. Mass Function 

Similar to the age estimation, the counting in between the 
CMD tracks was adopted to determine the distribution of 
masses. Nine evolutionary tracks ranging from O.1M to 7 
M were used to estimate the mass of pre-MS stars. For the 
MS objects we adopted the Padova model that best fits the 
observed colours, as described in Sect. 4.2. By this way, the 
sum of MS and pre-MS stars corresponds to the number of 
observed objects, which is given in Table 2 along with the 
respective observed mass. 

Instead of using a histogram to display the mass 
distribution, w e calculate the mass function given by 
Kroupa (2001) 



-(i+x) 



(6) 



Figure 3 (right panel) displays the observed mass dis- 
tribution based on the sum of both MS and pre-MS stars. 
By fitting the observed distribution, we obtained the slope 
of the mass function, represented by \ ( see ^ as ^ column of 
Table 2). 

About half of the sample has a mass distribution flatter 
than the initial cluster mass function (ICMF) verified in a 
large variety of clusters (e.g. |Elmegreen 2006| |Oey 2011] ). 



The ICMF slope ( X ~ 1.0) is slightly shallower than 

Salperter's IMF ( X ~ 1.35). 

The IMF suggested by |Kroupa (2001) assumes slopes 
X = 0.3 ± 0.5 for M < O.5M and Salperter's IMF for 
M > 1 M . We used the the observed mass and number 
of members to estimate the average stellar mass (m) and 
verified that most of them have fh > 1 M , probably due 
to a lack of low mass stars in our sample. This is related 
to the 2MASS detection limit that constrains the presence 
of faint sources in our sample. A sample incompleteness 
gives for our clusters large m, when compared to embedded 
clusters. Bonatto & Bica (|2010[ I20TT]) , for instance, uses fh 
— 0.6 M , which is lower than our results by a factor of 2, 
at least. 

Based on the IMF fitting of the observed distribution 
of mass (x slope), we corrected the incompleteness of our 
sample by synthetically deriving the number of faint stars 
that should be considered in the real membership. Table 
2 gives the total number of members Nx = N 5 S + N^ m /, 
where N 5 S is related to the observed objects and N^ m j is 
the number of lacking faint stars, estimated by integrating 
the IMF in the range of low-mass (below the limit of detec- 
tion). The corresponding total mass is also given in Table 
2, along with the observed mass. The correction on the 
sample completeness is needed to improve the estimation 
of the cluster parameters that are used in the comparison 
with other samples. 



5. Comparative Analysis 

We investigate possible correlations among the cluster char- 
acteristics by comparing the structural and fundamental 
parameters, respectively from Table 1 and Table 2. 

First, the correlations among mass, radius and number 
of members are compared to similar relations obtained by 
AdamseiaL](2006) and Pfalzner (2011) using cluster prop- 
erties from other works. In Sect. 5.4 the structural param- 
eters (obtained on basis of surface density) are compared 
among each other, and Sect. 5.5 presents the relations with 
age. 



5.1. Radius 

Figure 4a shows the cluster radius distribution as a function 
of number of members (the same as fig. 2 from I Adams et aLl 
20 06). Our sam ple is compared with 14 stellar clusters listed 
by |Carpenter| (2000), and 34 embedded clusters from lLL03l 
which have available size (radius in pc). 

In order to illustrate the different cluster size defini- 
tions discussed in Sect. 3.2, we also plot the cluster ra- 
dius (R t ) and the respective number of members (Ni) that 
were obtained by fPiskuno v et al.l for three of our clusters 
(NGC 2362, NGC 6604, and Stock 13). It can be noted that 
the relation R t x N± (crosses in Fig. 4a) is incompatible with 
the distribution of the other samples, possibly due to the 
fact that Kharche nko et al.l (2005) and citePiskunov 2007 
study different stellar groups that encompass our clusters 
(see Sect. 3.2). 

We verified that 52% of our sample follows the depen- 
dence R =0.1N 5 that is the same relation proposed by 
I Adams et al.l (2006) for the clusters listed by Carp enter] 
(2000). Two of our clusters, NGC 6604 and Trumpler 18, as 
well as three clusters from 1LL03I fNGC 2282, Gem 1, and 
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Table 2. Fundamental parameters. 



Cluster 


Age 


E(B-V) 


d 


(m-M) 


ffc 


M 


ob 




N obs 


M T 


N T 




X 




Myrs 


mag 


(PC) 


mag 


(%) 


M 




(stars) 


M 


(stars) 




Collinder 205 


5 +b - u 


0.56+0.07 


isootSX 


11 o + U.b 
ii -^-0.7 


14+ iY 

14 -10 


390 


± 


79 


174 ± 38 


876 


± 262 


2617 ± 


701 


1.17+0.33 


Hogg 10 


3 +2 -° 


0.23+0.05 


2200±I° ° 


n 7 +o.6 

± - L - ' -0.4 


99 + I8 
ZZ -15 


218 


± 


43 


88 ± 10 


318 


± 81 


465 ± 


149 


0.59+0.17 


Hogg 22 a 


3 +2 -° 


0.43+0.07 


1700 +550 


11 9 +0 - 6 




285 


± 


47 


97 ± 33 


376 


± 99 


429 ± 


220 


0.54+0.25 


Lynga 14 b 


3 +2 -° 


0.50+0.10 


950if ° 


Q Q + - 4 
y ' y -0.5 


e±S B 


127 


± 


30 


68 ± 16 


374 


± 111 


1374 ± 


430 


1.24+0.21 


Markarian 38 


5 5 +4 - 5 


0.10+0.07 




10 9 +0 ' 3 


21+ 11 


137 


± 


28 


57 ± 3 


190 


± 57 


251 ± 


107 


0.52+0.27 


NGC 2302° 


9 c + 2.5 
Z -°-2.0 


0.19+0.05 


1700l»g 


11.2+°; 3 6 


26+i? 


168 


± 


32 


70 ± 19 


238 


± 57 


345 ± 


111 


0.60+0.19 


NGC 2362 


40 + 6.0 


0.06+0.05 


1480±« 2 °° 


10.9+-1 


Q1+17 

CJ± -13 


332 


± 


55 


124 ± 25 


500 


± 126 


927 ± 


331 


0.92+0.25 


NGC 2367 


3 +2 - 5 


0.15+0.05 


2200±J°° 


11 7 +0 - 4 

L1 - 1 -0.2 


20+ 12 
zu_ 15 


155 


± 


29 


60 ± 6 


267 


± 64 


586 ± 


167 


1.03+0.21 


NGC 2645 d 


5 +5 -° 


0.28+0.08 


1800± 8 2 °° 


11 Q + 0-8 
ii -^-0.3 


27lf 8 


238 


± 


49 


104 ± 14 


498 


± 102 


1366 ± 


265 


1.10+0.20 


NGC 2659 




0.25+0.05 


2000± 2 ™ 


11 5 +0 - 5 




495 


± 


103 


215 ± 27 


857 


± 237 


1801 ± 


608 


0.89+0.22 


NGC 3572 


3.0±™ 


0.15+0.07 


1900±«S 


11 4 +0 - 4 

- L - L - -0.2 


28±1S 


371 


± 


72 


149 ± 22 


645 


± 175 


1422 ± 


485 


0.97+0.26 


NGC 3590 




0.35+0.07 


168011™ 


11 l + °- 6 


34+23 


175 


± 


36 


79 ± 20 


328 


± 90 


768 ± 


264 


0.95+0.27 


NGC 5606 e 


2 5 +2 - 5 


0.32+0.05 


2200+ 300 


n 7 +o.3 

J " L - ' -0.2 




262 


± 


49 


98 ± 11 


417 


± 81 


768 ± 


147 


0.86+0.13 


NGC 6178 f 


5 0+ 2 - 5 


0.15+0.08 


14 o +l50 


10.7+°- 3 


6+ 5 


220 


± 


52 


106 ± 7 


379 


± 284 


1430 ± 


1345 


1.43+1.38 


NGC 6604 


6 +4 -° 


0.67+0.07 


1600±S°S 


11 o+° :5 




279 


± 


52 


90 ± 7 


326 


± 95 


245 ± 


141 


0.26+0.21 


NGC 6613 


5 0+ 5 - 


0.45+0.07 


1550±»SS 


11 o +0 - 4 


26±iS 


345 


± 


73 


133 ± 9 


557 


± 113 


1087 ± 


185 


0.93+0.13 


Ruprecht 79 


5 +1 ° 


0.28+0.05 


2700± 4 3 ° 


-, 9 9 + 0.3 
1Z - Z -0.3 


40+ 20 

^ U -25 


459 


± 


86 


174 ± 8 


677 


± 119 


1156 ± 


154 


0.95+0.14 


Stock 13 


4 o+ 6 - 

^• U -3.0 


0.07+0.05 


2000±!|°° 


11 5 +0 - 5 


09+I8 


179 


± 


31 


65 ± 11 


308 


± 80 


686 ± 


236 


1.02+0.26 


Stock 16 


7.0+5"? 

— 6.5 


0.28+0.05 


2000^°° 


11 5+ - 5 
iJ -°-0.4 




254 


± 


66 


139 ± 13 


1043 


± 239 


4676 ± 


995 


1.50+0.18 


Trumpler 18 


5 0+ 5 - 


0.10+0.07 


2850± 4 3 °° 


-, 9 O + 0.3 
iZ -^-0.4 


26±« 


512 


± 


85 


164 ± 8 


572 


± 214 


334 ± 


543 


0.00+0.17 


Trumpler 28 


2 +3 -° 

Z ' U -1.5 


0.65+0.05 


1050lf ° 


10 i +0 - 5 

1U -0.5 


iol 25 


178 


± 


35 


73 ± 7 


233 


± 58 


268 ± 


86 


0.42+0.16 



Columns description: (1) Identification; (2) age; (3)colour excess; (4) Distance; (5) distance modulus; (6) fraction of stars showing K excess; (7, 8) 
observed mass and number of members; (9, 10) total mass and number of members; (11) mass function slope. 

Notes: Clusters showing a second peak on age distribution: (a) 12. 5+ 2 ^; (b) 12.0l 3 ;°; (c) 12. S^ 2 ^; (d) 12.0+ 3 ;°; (e) 11.5^1^ (f) 12. 5+ 2 ^. 



Gem 4) also follow this trend, but scaled u p by a factor 
1.7, which is the superior limit suggested by I Adams et al.l 
(2006) for the distribution of cluster size as a function of 
N. 

Coinciding with most of the ILL03I clusters, 33% of our 
sample are distributed along the lower curve in Fig. 4a that 
is similar to the |Carpenter| (2000) sample, but scaled down 
by a factor 1.7. 

Lynga 14 appears bellow this correlation, suggesting 
that its membership is larger than the expected number 
of members, when compared to other clusters having same 
size. This characteristic is also noted for some clusters from 
ILL03[ in particular S 106. 

We also evaluate for our sample the commonly used 
half-mass radius 7*1/2, which encloses half of the total mass 
of the cluster. Considering that we do not determine indi- 
vidual mass for each member of the clusters, their observed 
radial mass distribution cannot be established. In order to 
have an approximate estimation of 7*1/2, we adopted the 
integrated mass distribution M(r) given by Ada ms et al.l 
(2006): 



M(r/r Q ) 

M tot 



(r/r Q ) a 



(l + (r/r )) a 



(7) 



where r Q is the scale length that we assume to be the radius 
of the cluster. The validity of using Eq.7 for our sample was 
checked for Lynga 14, Collinder 205 and Hogg 10, for which 
we could estimate the observed M(r) and were used as test 
cases. The "virial" model, used by Ad ams et al.l (2006) in 
the simulations for N=100 members, has coefficients (a=3 
and p ~ 0.41) that seem to well reproduce the radial profile 
of the checked clusters. We verified that these simulations 
provide a ratio of ri/ 2 /r of about 60-63%. 

By adopting the relation 7*1/2 = 0.615 r G , we obtain the 
range of 0.69 - 2.5 pc for the half-mass radius estimated 



in our sample . We compared this result to the study that 
lAdams et aLl f2006) developed for NGC 1333, which model 
is most like the N=100 simulations with "cold" starting 
states (a=2 and p ~ 0.55) and gives 7*1/2 = 0.117 to 0.238 
for r o =0.3 to 0.4 pc. The relation ri/ 2 /r is about 40-60% 
in this case, which is probably due to the fact that our 
clusters are more evolved than the embedded phase. 

5.2. Mass 

Figure 4b shows the relations between mass and number of 
members. The LL03 sample has a dependence M = a N 10 , 
where a = 0.6 io*3> wn i cn means that they are distributed 
between the lines scaled (up and down) by a factor 2. Our 
objects also follow this dependence. 

It is interesting to note that the clusters NGC 2282, 
Gem 1, and Gem 4, from lLL03l are found above the superior 
limit of the distribution. The same occurs for our clusters 
Trumpler 18 and NGC 6604, which showed a similar trend 
in the R x N plot (Fig. 4a). 

The reason for a deviation of the expected relation is 
that these objects are different of most of the clusters. They 
have a smaller number of low-mass members, which is con- 
firmed by their flat IMF (x < 0.4). In the opposite sense, 
NGC 6178 and Stock 16 are found below the lower limit 
curve, being compatible with their steep IMF (x > 1.4). 

In order to discuss if these differences co uld be inter- 
preted as formation condition, as proposed by Ada ms et al.l 
(2006), or differently if they are related to a time sequence, 
as suggested by Pfalzner (j2011|) . we compare the relations 
between mass and radius. The thick line in Fig. 4c shows 
that our samp le has the same dependence obtained by fit- 
ting the QXOI data: M = 118 R 1 3 . 

Figure 4c also displays the data of the "leaky" mas- 
sive clusters studied by Pfalzner (2011), who proposed 
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R(pc) 



Fig. 4. Comparing our objects (filled circles) with embedded clusters (open squares) studied bv lLLQ3l and other samples. 
Representative error bars are used in a few data points, for illustration, (a): Radius of the cluster as a function of number 
of members. Our sample, as well as t he embedded clusters follow the same dependence R cx TV 5 (th ick line) found for the 
stellar clusters (asterisk s) studied by |Carpenter| (2000). Thin lines indicate the limits suggested by I Adams etHl (2006). 
The results obtained by Pisku nov et al.l (2007) for three objects of our sample (crosses) are used to illustrate differences 
on the cluster size definition, (b) : Mass of the cluster as a function of number of members showing the same dependence 
M cx N 1 for our sample and embedded clusters (thick line), with limits scaled by a factor ~ 2 (thin lines), (c): Mass- 
radius dependence showing a mean distribution of M ~ 118 R 1 3 (tick line) that spreads by a factor 2 (thin lines). The 
parameters of massive ("leak") exposed clusters (triangles) are also plotted, as well as the dependence M = 359R 17 
suggested by Pfalzner (|201ip for the embedded clusters (dashed line), (d): Volumetric density as a function of radius, 
showing the dependence p = 28 R _1 7 presented by our sample and ILL031 data (thick line), limited by thin lines that are 
scaled by a factor 2. Dashed lines represent the relations proposed by Pfalzner (2011]). 



that these objects would be the ending of a time se- 
quence. However, the mass-radius dependence M ~ 359 R 1 7 
(dashed line), which Pfalzner (2011) used to illustrate the 
suggested time sequence, is steeper than the distribution of 
LL03 clusters. 



and S 106 (ILL03|h which differences we interpret to be more 
likely due to different formation conditions. 



5.3. Volumetric density 



Instead of confirming a sequence that ends in the ex- 
posed massive clusters, our results show that clusters from 
Pfalzner (2011) having radius > 10 pc follow the same M 
vs. R relation presented by ILL03l data and also our sample. 
On the other hand, the massive clusters with R < 10 pc are 
found above the upper line in Fig. 4c, similar to Lynga 14 



By assuming a spherical distribution, we estimated the vol- 
umetric density, defined by p = 3M/(47ri? 3 ), which is plot- 
ted in Fig. 4d as a functi on of r adius. For both samples: our 
clusters and those from ILL031 we found a similar relation 
p = 28i?~ 1 ' 7 . This is quit e similar to the results obtained 



by Camargo et al. (2010)| that find p cx R 



-1.92 
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Pfalzner (|201ip used a di fferent dependence p ~ 
100R~ 13 to represent the ILLQ31 data, that were compared 
to the relation p ex R~ A obtained for the massive exposed 
clusters (dashed lines in Fig. 4 c). However, the distribution 
of our clusters, as well as LL03 data, does not agree with the 
mass vs. density relation proposed by Pfalzner (2009, l20TT]) . 
As suggested in Sect. 5.2, the large massive clusters seem 
to follow the same trend of our clusters, while those with 
R < 10 pc are scaled up by a factor larger than 2. 

5.4. Surface density 

Were analysed the structural parameters looking for cor- 
relations among the clusters themselves, and to verify a 
possible relation between the cluster and its environment. 
Hatched areas in Fig. 5 illustrate the trends that were found 
for most of the clusters. 

Excepting for Trumpler 28, the observed average den- 
sity of the cluster (< n* >) increases with the background 
density (crbg), as shown in Fig. 5a. This indicates that dense 
clusters are found in dense background fields. 

Considering that N increases with R (see Fig. 4a) it 
is expected an anti-correlation of < n* > vs. R, which is 
indeed noted in Fig. 5b. By consequence, cr^g also appears 
anti-correlated with R. 

Core density (<to) is another structural parameter that 
also diminishes when cluster radius increases, but a consid- 
erable dispersion is seen in Fig. 5c, where almost six clusters 
are out of the observed trend. 

In spite of the fact that R is not expected to be related 
to the cluster distance (D), Fig. 5d shows a trend of R 
increasing with D. This is not surprising, since the similar 
angular sizes of our objects lead to large linear size (given 
in parsec) for more distant clusters. 

We plot in Fig. 5e the anti-correlation between core ra- 
dius (r c ) and <7o that is observed for most of the clusters, 
excepting Lynga 14 and NGC 3590. The core parameters 
also showed some trends when compared with other struc- 
tural parameters, like the correlation between <jo and 5 C 
(density contrast); and r c /R (ratio of core-size to cluster 
radius) decreasing with the raise of cjq. 

The evaluation of structural parameters is unprece- 
dented for most of our objects, but our results are com- 
parable with those available in the literature for other 
clusters. We verified that the ranges of values obtained 
by us are compatible with several kind of young clusters 
like: NGC 6611, an embedded and dynamically evolved 
cluster dBonatto, Santos Jr. fe Bica 2006[); NGC 4755, a 
post-embedded cluster (Bonatt o et al. 20 06); the low-mass 
open clusters NG C 1931, vdB 80, BDSB 96 and Pismis 5 
(jBonatto & Bical 2009a); NGC 2239, a possible ordi- 
nary open cluster ( jBonatto Bica, 2009b); and the dis- 
solving clusters NGC 6823 ( |Bica, Bonatto L Dutra, 2008[ ), 
Collinder 197, vdB 92 (Bonatto fc Bic a, 12010]) and 
Trumpler 37 (Saurin, Bica & Bonatto. I2012p . 

5.5. Age 

It would be expected that do should decrease with age once 
the members of older clusters could have had time to dis- 
perse, diminishing its surface stellar density. However no 
clear trend is observed in Fig. 5f, which displays different 



symbols to indicate objects that have two options of age, 
but only one of them was adopted. 

The fraction of stars having large E(H-K), fx, is ex- 
pected to be related to age because of the colour excess 
in the K-band traces circumstellar matter, which is indica- 
tive of youth, as mentioned in Sect. 4.1 (ILada et al. 1996). 
Figure 5g shows the distribution of fx as a function of age, 
with a dotted line indicating the fx = 20% limit, above 
which objects younger than 5 Myr are supposed to be 
found. Most of our objects have fx > 20%, but three of 
the youngest clusters are below this limit, showing no cor- 
relation of age with fx- 

Probably this lack of correlation is due to the fact that 
some of our clusters have a mixing of populations, which 
means: part of the members is < 4 Myr, while others are 4 
- 10 Myr. 

On the other hand, the lack of correlation is also found 
for Stock 16 that has fx = 30% and mean age of 7.0^.1'® 
Myr. An explantion is a possible field-stars contamina- 
tion in the number of objects with K-band excess in the 
colour-colour diagram, mainly those having 0.1 < (J-H) 
< 0.3 mag. In the colour-magnitude diagram, these stars 
are counted in the range of ages older than 20 Myr, but 
should not be considered. The large error bars on both, age 
and fx estimation, give us only a qualitative analysis of the 
relation between these parameters. 

Age is also expected to be related with other parameters 
such as colour excess. Since E(B-V) is indicative of visual 
extinction, it can be used to infer how embedded the cluster 
is, and by consequence to verify the youth of the cluster. 
In fact, Fig. 5h shows three of the youngest objects (< 5 
Myr) appearing above the line representing E(B-V) =0.3 
mag, which means Ay > 1 mag. However, five other young 
clusters are below this limit, which indicates that they are 
not deeply embedded. For this reason, we concluded that 
there is no correlation between age and E(B-V) for our 
sample. 



6. Summary of the results and conclusions 

We determined the structural parameters as a function of 
superficial density and radial distribution profile of a sam- 
ple of 21 clusters, selected on basis of their youth and in- 
termediate distance. A statistical procedure using colour- 
magnitude criteria provided a double-checked decontami- 
nation of field-stars. The remaining stars were considered 
members and their observed colours were used to more ac- 
curately determine the visual extinction affecting the clus- 
ter. The unreddened colours were compared with theoreti- 
cal isochrones in the CMD aiming to confirm distance and 
age. The same was done to determine individual masses, 
based on the cluster members position compared to evolu- 
tionary tracks. 

In principle, centrally concentrated clusters should be 
the youngest ones since they would not had enough time to 
disperse. In fact, most of our clusters show this character- 
istic. However, the constrained range of ages in our sample 
impeaches us to be more conclusive about differences or 
similarities on the evolution of the studied clusters. 

We conclude that all the 21 studied clusters are very 
similar, probably due to the selection criteria choosing re- 
stricted ranges of size, distance and age. By consequence, 
there is no large variation on number of members, radius 
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Fig. 5. Comparing structural parameters: (5a) background density vs. average density. Cluster radius compared to (5b) 
average density; (5c) core density; and (5d) distance. (5e) core radius vs. core density. Age correlations with: (5f) core 
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symbols correspond to the first choice, listed in Table 2). 



and mass of the clusters. On the other hand, the galac- 
tic distribution of the objects causes differences among the 
clusters environments. 

When compared with other young clusters (|LL03| 
|Carpenter| 2000), our sample follows the same trends, but 
has larger radius and lower volumetric density. This means 
a less concentrated distribution of members that may be 
related to the expected spatial dispersion, when the cluster 
gets older. 

The distinction between star clusters and associations 
has been discusse d in se veral works. Following the defini- 
tions presented bv lLL03[ our clusters are classified as stellar 
groups because they have more than 35 physically related 
stars and their mass density exceeds 1 M /pc -3 . Since our 
objects are optically visible, they cannot be considered em- 
bedded clusters. Even considering the large error bars on 
the age estimation, the mean age of the clusters (~ 5 Myr) 
is a clear indication that our sample is formed by young 
stellar associations. 

We also checked for our sample the relation 
between age and crossing time (r cr ), following 
|Gieles fc Portegies Zwart (2Q11)[ for instance. They 
proposed a criterion in which bound systems (open star 
clusters) would have age/r cr > 1 and unbound systems 
would have age/r cr < 1. For this test, we adopted the 
crossing time defined by: r cr = where R is the cluster 
radius and the stellar velocity dispersion is given by: 

ay = \J according to Saurin, Bica & Bonatto (|2012p . 

We verified that none of our clusters has stars with age 
that exceeds the crossing time. Therefore, they probably 
would evolve like stellar associations. If the error bars on 



the age/r cr calculations are considered, only Lynga 14 could 
evolve as an open cluster. This result is in agreement with 
the suggestion by LL03 and Pfalzner (2009), for instance, 
that few star clusters are expected to be bound. However, it 
must be kept in mind that these are qualitative conclusions, 
due to the large uncertainties on age estimation. 

Our sample fulfils the caveat between the samples of 
embedded and exposed clusters studied by Pfalzner (|201ip . 
However, our results do not confirm the mass-radius or 
density-radius dependences suggested by her. 

In fact, we verified that massive clusters from Pfalzner 
(|201ip are distributed in two groups. Those with large sizes, 
lower masses, and intermediate ages (4-10 Myr) follow the 
same trends shown by embedded clusters (LL03), as well as 
our sample. The younger massive clusters (< 4 Myr) stud- 
ied by Pfalzner (|201ip have lower sizes and higher masses, 
appearing out of the correlations shown in Fig. 4, as well 
as Lynga 14 (this work) and S 106 (LL03). As proposed by 
Ada ms et al.l (2006), the differences among these clusters 
could be interpreted as formation condition. 

An interesting perspective of the present work is to in- 
crease the studied sample by including other clusters hav- 
ing larger radius (4 - 10 pc), which could complet the gap 
between our clusters and the sample of massive clusters. 
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Appendix A: Plots of the entire sample 

All plots used in the analysis of structural and fundamental 
parameters are displayed in this Appendix. Figures Al and 
A2 show the stellar surface-density maps and the distribu- 
tion of the stellar density as a function of radius (same as 
Fig. 1). 

Figures A3 to A7 present colour-colour and colour- 
magnitude diagrams in the left panel (same as Fig. 2). The 
centre and right panels show respectively the histogram of 
age and the mass distribution (same as Fig. 3). 
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Fig. A.l. Left: Stellar surface-density map (a (stars/arcmin 2 )) obtained for the region of 30 arcmin around the clusters. 
The comparison field-stars area is indicated by dashed lines, while the full line indicates the cluster area. Centre: a zoom 
of the <T map indicating by crosses the position of objects with membership probability P>70%. Right: The distribution 
ht the stellar density as a function of radius. The best fitting of observed radial density profile, indicated by the full line, 
was obtained by using the model from King (1962). A dashed line indicates the background density ((Tb g ). 
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Fig. A. 2. Same as Fig. Al 




Fig. A. 3. Left panels: Colour-colour and colour-magnitude diagrams. The MS and the ZAMS are indicated by full lines, 
while the locus of giant stars is represented by a dotted line. Reddening vectors from Rieke & Lebofsky (1985) are 
shown by dot-dashed lines. The isochrones and evolutionary pre-MS tracks from Siess et al. (2000). Cluster members are 
indicated by open diamonds (P50) and dots (P?). Centre: Age distribution of clusters members (thick line) showing the 
contribution of P50 (dotted line) and P? (dashed line) objects. Right: Observed mass distribution indicated by crosses 
with error bars. The thick line represents the mass function 0(m) fitting.) 
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. A. 4. The same as Fig A. 3. 
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. A. 5. The same as Fig A. 3. 
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. A. 6. The same as Fig A. 3. 
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. A. 7. The same as Fig A. 3. 



